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Abstract 



We investigate Fano-type anomalous transmission of energy of plane waves across lossless slab scat- 
terers with periodic structure in the presence of non-robust guided modes. Our approach is based on 
rigorous analytic perturbation of the scattering problem near a guided mode and applies to very general 
structures, continuous and discrete. 

1 Introduction 

Resonances with sharp peaks and dips have been observed in the graphs of the frequency dependence of 
the transmission of energy of plane waves through periodic slab structures. GlElEEOlO]]. These anomalies 
have been connected to isolated real points on the complex dispersion relation D(k, u) = that relates the 
frequency u of a sourceless field to its Bloch wave vector k along the slab (H [TOl . The real point of the 
dispersion relation corresponds to a true guided mode in the slab, that is, one that decays exponentially away 
from the slab. Because the point is isolated, the mode is non-robust with respect to perturbations of k, lo, or 
the physical or geometric parameters of the structure. This feature gives rise to sharp transmission anomalies 
near the frequency of the guided mode, as shown in Figs. 1 and 2. 

The anomalies observed resemble those investigated by U. Fano [4] in the context of quantum mechanics, 
in which an embedded eigenvalue is "dissolved" into the continuous spectrum upon perturbation of the 
system, resulting in resonant dynamical behavior. The shape that he derived for the resonance involves two 
parameters and turns out to be a special case of a more general shape derived rigorously for the graph of 
energy transmission through photonic crystal slabs with a symmetry assumption on the slab geometry ifTTTl . 
Typically, symmetry gives rise to standing guided waves (k = 0), as discussed in |E1 U2J, and only standing 
waves were analyzed in fH). m the present work, we compute anomalies near traveling guided modes 
(k ^ 0). The frequency of a non-robust guided mode is an embedded eigenvalue for a fixed value of n. 
The dissolution of the eigenvalue into the continuous spectrum corresponds to the destruction of the guided 
mode when n is perturbed. We demonstrate the results through two examples, the scattering of polarized 
EM fields by a lossless dielectric structure and the scattering of waves in a two-dimensional lattice by a 
one-dimensional periodic lattice attached along a line. 
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2 Analytic connection of scattering states and guided modes 



In ifTTIl . we reformulate the Maxwell equations in a photonic crystal slab (Fig. 2) at fixed frequency oj and 
Bloch wavenumber k (the latter necessarily parallel to the slab) as a boundary-integral equation (see HQ, 
for example). The equation 

A(k, oj)ij) = (j) (1) 

relates two fields and ip defined on the interface between the scatterer material (taken to be homogeneous 
and isotropic) and the ambient medium. The field tp on the interface contains information that is necessary 
and sufficient for the calculation of the EM field in the whole space; the calculation consists of the evaluation 
of an integral. Similarly, the field 4> contains information for the calculation of the incident field. Solving the 
above equation for i/> is equivalent to calculating the field in the photonic crystal and in the ambient medium 
established as a result of the source field cj>. We call ?/> the total field to distinguish it from the scattered field, 
which is the difference between the total and incident fields in the ambient medium; clearly, a scattered field 
is required to satisfy a radiation condition. The integral operator A is of second-kind Fredholm type with 
index zero, thus, it either has a bounded inverse or it has a nullspace of finite dimension and a range of the 
same codimension. 

The integral operator A is parametrized by oj and k. When either of these takes a nonreal value, a field ip, 
bounded on the interface, may produce an unbounded and hence unphysical field in space due to growth of 
the exponential e lKX or e lujt . Such fields are not only mathematically but also physically useful and form the 
foundation of the so-called leaky modes. A (generalized) guided mode is a nonzero solution of the sourceless 
equation 

A(k,u)i/) = 0. (2) 

The pairs (k, oj) for which such a solution exists satisfy a dispersion relation D(k, oj) = 0. For pairs («, oj) 
such that k is real and D(k, oj) = 0, we must have Imw < lHJEEQl. In the case that Imw = 0, a solution to 
(|2]) represents an EM field that decays exponentially with distance from the scatterer and therefore represents 
a true guided mode. We shall work with a simple eigenvalue branch £(k,oj) of A(k,oj), for which we 
necessarily have that £(k, oj) = implies D(k, oj) = 0. 

The framework we have described explicitly for the Maxwell equations arises for very general dielectric 
and metal structures (not necessarily with homogeneous components) and for other continuous and discrete 
problems of scattering by a slab; the methods of arriving at the formulation vary from problem to problem. 
For waves in an n-dimensional uniform lattice scattered by a (n — 1) -dimensional periodic lattice (Fig. 1), 
the problem is reduced to a finite-dimensional one, and A is a matrix. 

For the purpose of investigating anomalous scattering behavior near a guided mode, we let (kq, ojq) be a 
real pair that satisfies £(kq, ojq) = and use as an incident source field the plane wave 

(j) inc (x, z) = £(k, oj)e iKX e ir > z (incident from the left), (3) 

where x is the multi-variable in the directions parallel to the slab and z is the single variable perpendicular 
to the slab (x and z may be discrete or continuous). The number 77 depends on k, oj, and the structure itself. 
We shall work in the regime of only one propagating diffractive order, that is, one Fourier harmonic in the 
periodic variable x that carries energy in the z-direction. In this case, the reflected and transmitted fields (at 
far field) are simple: 

<p iefl (x, z) = a(K, oj)e iKX e~ iriz (reflected on the left), (4) 
<j)tmns{x,z) = b(K,oj)e lKX e tT > z (ttansmitted on the right). (5) 
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It can be shown that the coefficients a(«, u) and lo) are analytic functions of (k, lo) near («o, loq) EMU]. 
At pairs (k,u;) that satisfy £(k,lo) = 0, the the right-hand side of ([I]) is zero, and the solution is a guided 
mode. Thus we obtain an analytic connection between scattering states and generalized guided modes near 
the true guided mode at (kq, loq). 



3 Analysis of transmission near a guided mode frequency 

The analysis of the resonant transmission shape near a guided mode is based on the analysis of the three 
analytic functions £(k,lo), a(n,uo), and u) near the real parameters (kq,loo) of a guided mode. The 
mode condition 1(ko,luo) = 0, together with that of conservation of energy, namely \£\ 2 = \a\ 2 + |6| 2 for 
real values of k and lo, implies that (kq , loq) is a root of the three analytic functions £, a, and b simultaneously: 

£(k ,lo ) = 0, a(K ,wo)=0, b(K ,LJo)=0, at (kq, loq) G M 2 . (6) 

The perturbation analysis of £, a, and b near this common real root relies principally on the following two 
conditions, which we have already mentioned: 

\£(k,lo)\ 2 = \ci(k,oj)\ 2 + \b(K,co)\ 2 for real k and lo, (7) 
if £{k, lo) = for k e M, then Im uj < 0. (8) 

The following conditions are assumed because they hold generically: 

d£ da 3b 

— (k ,u ) / 0, — (k ,u;o) / 0, t^(ko,wo)/0. (9) 

The WeierstraB preparation theorem for analytic functions of several variables [5 ] then provides the following 
forms (with w = u — ujq and k = k — kq): 

£{k, lo) = e idl [w + I x k + £ 2 ~k 2 + 0{k 3 )] [1 + 0(\R\ + \w\)\ , (10) 
a(K,uj) = e ie2 [w + n« + r 2 k 2 + 0(k 3 )] [r + 0(\k\ + \w\)} , (11) 
b(K,u) =e i03 [w + hk + t 2 k 2 + 0(k 3 )] [t + 0(\k\ + M)] . (12) 

The function £ has been normalized so that the constant in the second factor is equal to 1 . All factors in these 
forms are analytic, and the numbers ro and to are real and positive. 
Consequences of ([8]) are that £\ is real and Im £ 2 > 0: 

r > 0, t > 0, 4 G M, Im£ 2 > 0. (13) 

Using the forms ( 10|11|12 1, we can examine various terms in the expansion of the equation of conservation 
of energy ((7]) for real values of k and lo. Computation of |£| 2 yields 

££ = [ w 2 + £ 2 k 2 + 2£ x wk + 2Re £ 2 wk 2 + 24Re £ 2 k 3 + 

+ (2^Re4 + \i 2 \ 2 )fi 4 + ...] [1 + 0{\k\ + \w\)\ , (14) 

and computation of \a\ 2 and \b\ 2 yield analogous expressions. In the case that £\ / 0, three terms of the 
expansion of |£| 2 = \a\ 2 + \b\ 2 are determined in terms of the coefficients expressed explicitly in the forms 
( |1Q|1 1|1 2] ) . In the case that £\ = 0, some of these terms become trivial and others become relevant. We 
analyze these two cases separately. 
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3.1 Case 1:1^0 



The three terms in the expansion of (|7]) that give information about the coefficients in the forms ( 10|11|12 1 
are 



w term: 1 = r$ + t, 



R 2 term: 



Hjhl 2 + *g|ti| 2 , 



?! =rgReri + tgReti. 



ctk term: 

Because of the convexity of the function x x 2 , the ro 2 -term and the ru/i-term imply that 

t\ < r 2 (Re n ) 2 + t 2 (Retx) 2 , 
with equality if and only if Re r\ = Re t\, and the K 2 -term is expanded to 

t\ = r 2 (Ren) 2 + t 2 (Retx) 2 + r 2 (Imn) 2 + t 2 (Im^) 2 . 



From conditions ( fT6| ) and ( [17] ) together, we infer that 

Imri = Imti = and t\ = t\ = l\. 
We thus obtain useful expressions for the zero-sets of I, a, and b near the guided mode pair (kq, loo) 



£(k, uj) 
a(n, uj) 
6(k, uj) 









u = uo — £±(k — Ko) 
a; = wo — — ^o) 
CJ = ojq — — kq) 



1 2 (k 
i 2 (« - Ko)" 



«0j 
Ko) 



(n G E) 



(15) 

(16) 

(17) 

(18) 

(19) 
(20) 
(21) 



We observe that, if all coefficients r n and t n vanish, then the curve a(/t, uj) = (resp. b{n, u) = 0), for 
real values of k, describes real frequencies to for which transmission T is equal to 100% (resp. 0%). We 
do know that l\ is real, so that these curves lie at most a distance 0((k — ko) 2 ) from the real w-axis, and 
one can deduce that the real parts of the frequencies on these curves correspond to peaks and dips of the 
transmission, which do not necessarily reach exactly 100% or 0%. 

These results allow us to make two important observations about the shape of the transmission resonance 
as a function of frequency. Both of these are illustrated on the graphs in Fig. 1 . 

1. For k ^ Ko, either both the peak and dip of T as a function of uj lie to the left of the frequency ujq or 
both lie to the right of ojq. On which side they lie depends on the sign of l\ and k — Kq. 

2. The order in which the peak and dip in T occur on the real w-axis is the same for k < kq as it is for 
k > kq (assuming T2 / £2)- This is because the coefficients of the linear terms in ((20]) and ( |2T] ) are 
equal. If T2 < £2, then the peak comes to the right of the dip, and if t<i < ri, then the peak comes to 
the left of the dip. 

The transmission coefficient (square root of transmitted energy) as a function of real k and uj is 



T(k,u) 



b(n, uj) 



£(k, uj) 



to 



\w + l\k + t 2 k + • • • 

\zu + £lK + £ 2 k 2 -\ 



|1 + 771 w + 772K + 



to 



W + Ilk + t2K 



W + t\K + I2K 

The coefficients have the following significance: 



(1 + Re?7i w + Re 772 «) near(Ko,^o). (22) 
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1 . l\ is the rate at which the anomaly in u moves past the guided-mode frequency ujq as a function of k. 

2. lim T(k, uj) = to- m particular, T(k, uj) is continuous at the guided mode pair (/to, ujq)- 

3. J^T(kq,uj) |o, = toRe This relates Re 771 to the slope of the non-resonant transmission at k = 

' t 2 -h 



4. £T(«,uto)| w = toRe(a^ + f72 

5. The coefficients £2 and £2 control the spreading of the peak and dip of the anomaly as k is perturbed 
from kq. 
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Fig. 1. Plane waves in the uniform 2D lattice are incident upon the coupled periodic ID lattice from the left and are 
transmitted to the right. The square root of the percentage of transmitted energy as a function of frequency uj is shown 
for several values of the Bloch wave vector k near the parameters (kq, ujq) — (0.06167, 0.979167) of a guided mode. 
The exact calculation and the theoretical formula are practically identical. 



3.2 Case 2: 







This case occurs when the dispersion relation £(k, uj) = 0, solved for uj, is symmetric in the variable k about 
no. In particular, this occurs at resonant frequencies for ko = if the structure has symmetry about a line or 
plane perpendicular to the plane of the scatterer OHHOlOl, as is the case in Fig. 2. 

From the second relation in §15\ , we observe that £\ = implies r\ = and t\ = also, and we find that 
the following three terms in the expansion give information about the coefficients: 



w 2 term: 1 = r\ + t{), 

wk 2 term: Re £ 2 = r 2 . Re r 2 + t 2 Re t 2 , 

K 4 term: \£ 2 \ 2 = r 2 \r 2 \ 2 + t 2 \t 2 \ 2 . 



(23) 



Because £\ = 0, the expansions ( |19|20|21 1 reduce to 





uj) 


= 4= 


=> UJ 


= w 


-£ 2 (k- 
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2 - 


a(n 


uj) 


= <= 


=> UJ 


= w 


- r 2 {n - 


«0 


) 2 - 


b(n 


uj) 
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- t 2 (K - 
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2 - 



(24) 
(25) 
(26) 
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The second of the relations ( [23] > tells us that the curve £(k, uj) lies between the curves a(n, uj) = and 
uj) = 0. In particular, if r 2 and t 2 are real and £ 2 is imaginary, then the transmission peak and dip move 
in opposite directions away from the guided mode frequency ujq as k is perturbed from kq. 

In iTTn . we demonstrate that the following approximate expression for the transmission anomaly near 
(kq, ujq) matches numerical simulations for the scattering of ^-polarized fields by a periodic lossless dielec- 
tric slab: 

tl\w + t 2 k 2 \ 2 {l + rjw) 2 



T {K,UJ) ~ r V + r 2 R2|2 + t 2| CT + ^ 2 |2 (1+r?ro) 
Here, the coefficients have the following significance: 



near (kq, ujq)- 



(27) 



1. lim T(kq,uj) = to and lim T(k,ujq) = to 

UJ — >UJQ K — >KQ 



2. 77 



1 dT 

t r 2 Ouj 



uj - 



3. The coefficients t 2 and r 2 control the spreading of the peak and dip of the transmission anomaly as k 
is perturbed from kq. 




Real frequency (co) 

Fig. 2. Transmission of E-polarized EM fields through a periodic array of infinitely tall lossless rods (e = 12, fi = 1) 
as a function of the reduced frequency uj, for various of values of the wave vector n near the parameters (0, 0.669) of 
a guided mode. The exact calculation and the theoretical formula are practically identical. 



3.3 Relation to the original Fano shape 

One of the classic examples of anomalous scattering behavior is observed in the excitation of the noble gases 
near characteristic energies of the atom (the Auger states [9]). The anomalies exhibit a peak and a dip and 
are called "Fano resonances" after the work of U. Fano H, in which he derived a formula for this shape: 

a = const. i - j 2 , f= — , (28) 

where the parameters q and F control the locations of the peak and dip and / is the deviation of the frequency 
from resonance normalized to a characteristic width F. 
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There is a connection between formula ( |27| ) and the Fano shape ( |28| ) that can be expressed in concrete 
terms. Namely, the six real parameters in (27 1 (to, fj, and the real and imaginary parts of ti and r-i) can be 
reduced to two if certain conditions are satisfied, resulting in the Fano shape. The conditions are (in addition 

to ti = 0) 

1. T2 and t2 are real (meaning that the extremal values of T are and 1), 

2. 7] = (the background transmission is flat), 

3. rh-2 + tlt2 = (for small real R, the dispersion relation given by d24b is purely imaginary). 



The resonance ( |27] >, as a function of frequency u, now reduces to the Fano shape ( |28| > with 

r = 2«V(^o) 2 + (tto) 2 and 9 = t/yVo) 2 + («o) 2 • (29) 
The first condition is satisfied by the example in Fig. 2, while the second and third are not. 

4 Phase anomalies and resonant amplitude enhancement 



The expansions ( 19|20|21 1 can be utilized also to yield a formula for the shape of the argument of the complex 



transmission coefficient near the guided mode parameters (kq, u)q), 

phase of transmitted field = arg b(n,uj) — arg £(k, uj) , (30) 

which exhibits sharp spike as a function of u> for values of k close to kq. This quantity is closely related to 
the effective density of states associated with transmission of incident waves through a periodic structure HI . 

A phenomenon that is directly associated to the interaction of plane waves with guided modes is the res- 
onant amplitude enhancement of the scattered field inside the structure. The scattering problem is nearly 
singular for small perturbations, so the solution in the scatterer is large compared to the incident field and 
resembles a guided mode. We have established in ifTTTl a leading-order result for the behavior of this en- 
hancement as a function of k: 

const 

amplitude enhancement ~ — . — . — (k — ► 0). (31) 

| ft | 

One of the essential ingredients in the analysis is the fact that there is no resonant enhancement at k = kq. 
This fact is tantamount to the fact that the scattering problem has a solution at the parameters of the guided 
mode, although it is not unique [2]; in other words, the source field has no "resonant component" at normal 
incidence. 
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